Exact integral equation for the renormalized Fermi surface 
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The true Fermi surface of a fermionic many-body system can be viewed as a fixed point manifold 
of the renormalization group (RG). Within the framework of the exact functional RG we show that 
the fixed point condition implies an exact integral equation for the counterterm which is needed 
for a self-consistent calculation of the Fermi surface. In the simplest approximation, our integral 
equation reduces to the self-consistent Hartree-Fock equation for the counterterm. 
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In his authoritative book on interacting Fermi systems 
Nozieres wrote fourty years ago Q: "In practice, we shall 
never try to calculate the Fermi surface, which is much 
too difficult." What is the reason for this difficulty? For- 
mally, the Fermi surface of an interacting Fermi system 
is denned as the set of all wavevectors satisfying 0] 



£kj, - M + S(k F ,iO) = 



(1) 



where £k is the energy dispersion in the absence of inter- 
actions, [i is the chemical potential, and S(k, to) is the ex- 
act self-energy of the interacting system |jj . For simplic- 
ity we assume an infinite and spin-rotationally invariant 
system at zero temperature, so that S(k, oj) is indepen- 
dent of the spin. Unfortunately, the function £(kp,i0) 
in Eq. (Q) is not known a priori, so that the calcula- 
tion of the true Fermi surface requires the solution of the 
many-body problem. 

For weak interactions, one might try to determine 
the Fermi surface perturbatively by simply calculating 
E(k, zO) in powers of the interaction and substituting 
the result into Eq. (|l|). However, in general the per- 
turbation series contains anomalous terms B with un- 
physical singularities, which are generated because the 
ground state of the non-interacting system evolves into 
an excited state of the interacting system when the in- 
teraction is adiabatically switched on. As discussed by 
Nozieres Q, this artificial level crossing can be avoided 
by introducing counterterms which are determined by the 
requirement that the Fermi surface remains fixed as the 
interaction is adiabatically switched on. This intuitive 
idea can be implemented perturbatively as follows |) : 
Suppose we would like to know the true Fermi surface 
of a system with Hamiltonian H = Hq + Hi, where Hi 
describes some general two-body interaction and the non- 
interacting part is given by 



Ho 



£ kcL c k CT 

k,o- 



(2) 



Here Ck CT are the usual annihilation operators of fcrmions 
with momentum k and spin a. An expansion in powers 
of Hi leads to Feynman diagrams where vertices corre- 
sponding to Hi are connected by propagators Go(k, uj) — 




FIG. 1: Decomposition k = nkF(n) + p of a wavevector k 
into a component k^ = nkp(n) on the Fermi surface and 
a component p in the direction of the local Fermi veloc- 
ity vf- The thick solid line is a part of the Fermi surface. 
This construction defines kf and n as a function of k. Here 
vf = Vkek|k=kj7 is defined in terms of the gradient of the 
free dispersion at the true Fermi surface, so that vf is not 
necessarily perpendicular to the Fermi surface. 



[u> — ek + fi] • These are singular for u> = and 
£k = /A which is not the true Fermi surface defined in 
Eq. (|l|). If the perturbative expansion is truncated at 
a finite order, this leads to the unphysical divergencies 
mentioned above |Q. To avoid these, we add the coun- 
terterm X^ko- ^(k-Fj *0)c k(T Ck(T to Ho and subtract it again 
from Hi , writing H = Hq + H' x , with 



H'o 



[ek + S(k F ,iO)]cj c(T c k , 



(3) 



and H[ = Hi — J2ka ^(^f, iO)cj C(T Ckcr. Here kp is the 
wavevector closest to k lying on the Fermi surface, see 
Fig. |l|. Using Eq. ([!]), the corresponding free propaga- 
tor can then be written as Gg(k, u>) = [to — £k + £k F ] _1 , 
which by construction is singular on the true Fermi sur- 
face. But how do we find the counterterm necessary to 
calculate £k F ? Following the usual strategy adopted in 
field theory || , we may expand the irreducible self-energy 
associated with the modified interaction H[ perturba- 
tively in powers of H[ and require that, order by order in 



2 



perturbation theory, the corrections vanish when we set 
u> = and k = kp. This renormalized perturbation the- 
ory leads to complicated integral equations for T,(k F ,iO), 
which must be solved numerically 0. 

A few years ago Anderson || critically discussed a sim- 
plified version of the renormalized perturbation theory 
outlined above, where £(kp,zO) is replaced by a con- 
stant Sfi. He correctly pointed out that in general it is 
not allowed to ignore the momentum-dependence of the 
counterterms, and argued that in two dimensions the ef- 
fective interaction between electrons with opposite spin 
and wavevectors on the true Fermi surface depends in a 
subtle way on the boundary conditions which cannot be 
adequately taken into account perturbatively. To gain a 
better understanding of this problem, it should be use- 
ful to have an algorithm for calculating the Fermi sur- 
face which does not rely on the perturbative expansion of 
counterterms in powers of the interaction. In this work 
we show that such an algorithm follows in a straight- 
forward way from the exact functional renormalization 
group (RG) approach described in Ref. || . 

In the past decade several authors have used Wilso- 
nian RG methods to study interacting Fermi systems 
§, |0[ @ [jj> @ II: |H> usin g different versions 
of the RG. In particular, the RG transformations used in 
Refs. |l3[ [l5| |1| exclusively focus on the mode elimi- 
nation step. Although such a procedure is sufficient if one 
considers only the one- loop flow of marginal couplings, in 
general a complete Wilsonian RG transformation consists 
not only of the mode elimination, but includes also the 
rescaling of momenta, frequencies and fields p7| . While 
the field rescaling (i.e. the wavefunction renormalization) 
becomes only important beyond the one-loop approxima- 
tion, the RG flow of all relevant and irrelevant couplings 
is determined in an essential way by the rescaling step 
even at the one- loop level. We emphasize that the rescal- 
ing of momenta and frequencies is more than a trivial 
mathematical change of variables - it is crucial to detect 
possible fixed points of the RG and to calculate critical 
exponents |[7|. Moreover, as will be shown below, the 
length and time rescaling is very helpful for a derivation 
of the self-consistent Hartree-Fock approximation within 
the framework of the exact RG. The importance of rescal- 
ing in the presence of a Fermi surface has been empha- 
sized by Polchinski fnj and by Shankar f|. In Ref.|] 
we have shown how the exact functional RG approach de- 
veloped previously Jl3[ [l5| can be modified to include 
the rescaling step. Here, we shall show that the rescal- 
ing is crucial to obtain a non-perturbative definition of 
the Fermi surface as a fixed point manifold of the RG. 
We obtain an explicit algorithm for calculating the Fermi 
surface which does not rely on the iterative procedure of 
fixing the counterterms order by order in perturbation 
theory. The fixed point property of the Fermi surface 
has also been emphasized by Ferraz who recently 
discussed the Fermi surface renormalization in a special 



two-dimensional system using the field theoretical RG. 

The exact functional RG can be formulated in 
terms of an infinite hierarchy of coupled differen- 
tial equations for the irreducible 2n-point vertices 
rf n \K' l7 ... ,K' n ;K n ,... ,K 1 ), where K = (a,k,iu) is 
a collective label for spin projection c, momentum k, 
Matsubara frequency iui. Here £ is an infrared cutoff 
with units of energy which regularizes the singularity 
of the free propagator. To scale wavevectors toward 
the Fermi surface, it is convenient to perform a non- 
linear coordinate transformation in momentum space, 
k = nfci?(n) + Vi?£<7/|vi?|, and eliminate k in favor 
of the dimensionless variable q and the unit vector n. 
Here kp(ii) is the length of kp parameterized by n, 
and Vf = Vk£k|k F is the Fermi velocity of the non- 
interacting system at the true hp, see Fig. |l|. The 
corresponding unit vector is denoted by vj? = vp/|vp|. 
Geometrically, q = vp-p/^ = vp-(k — nfcp(n))/£ mea- 
sures the distance of a given k from the Fermi surface 
in units of the infrared cutoff. We also define rescaled 
frequencies e = ana - label the degrees of freedom by 
Q = (ex, n, q, ie) instead of K. To implement the scaling 
toward the Fermi surface [ fi"i"| we consider the RG flow of 
the following rescaled vertices M 



\{2n) 



(Q'i, 
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,Qi) 



1 £71 



1/2 



,K' n ;K n ,... ,K X ). (4) 

Here t = — ln(£/£o) is a logarithmic flow parameter, 
v ° = / (2 < f) p ( K £k — e k F ) is the density of states (per 
spin projection) of the non-interacting system at the true 
Fermi surface, and 



9S ? (kF ,uj + iO) 



(5) 



is the wavefunction renormalization factor, where 
S^(k, icu) is the irreducible self-energy of the system with 
infrared cutoff £ = £o e ~*- To determine the true Fermi 
surface self-consistently it turns out to be useful to sub- 
tract the counterterm E(kp, iO) from the irreducible two- 
point vcrtexfl, defining 



fi 2) (Q) = ^frf\K) = [E e (J0 - E(k F ,io)] 



zK(2) 



z? 



(6) 



The subtracted two-point vertex satisfies the exact flow 
equation (9) 



d t r[ 2) (Q) 



{l-r$-qd q -ed e )tf\Q) 



G*(g')fi 4) (Q,Q' ; g',Q) , (7) 



Q' 
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where 77" = — <9(lnZ" is the flowing anomalous dimen- 
sion (which vanishes for large t if the system is a Fermi 
liquid). The integration measure in Eq. (Q) is defined by 

= E/^/^(W^ - (8) 



■>D 



2 TV 



where dS^ is a surface element and Srj is the surface area 
of the unit sphere in D dimensions, and J(n,q) is the 
Jacobian associated with the transformation k — > (n, q) 
divided by v^ 2 . The function G t (Q) in Eq. (0) is defined 
by 



Gt(Q) 



s(n Q - 1) 



Z? [ie-&(qj\ +f| 2) (Q) 



(9) 



where £f(q) = (e kf -_ek F )/£, with k t = nfc F (n) + 
v F £oe~*g/l v -H- Here (Iq is some function that mea- 
sures the distance from the Fermi surface, for example 
Qq = |£ t n (<z)| ~ \q\- The flow equation of the rescaled 
irreducible four-point vertex f ^ (Q[, Q' 2 ; Q2, Qi) is not 
explicitly needed in this work; it can be found in Ref.j9). 

The shape of the Fermi surface of the interacting sys- 
tem is determined by the RG flow of the couplings 



flf =ffV,n l(? = CMe^0) 



(10) 



From Eq. (0) we see that these couplings satisfy the exact 
flow equation 



(11) 



with 



G t (Q')f, (4) (Qo,Q';Q',Qo) , (12) 



Q' 



where Qo = (tr, n, q = 0,ie = iO). Note that in dimen- 
sions D > 1 there are infinitely many couplings flf, la- 
belled by the unit vector n. Obviously, each jlf with 
n*^ = limt_>oo 17" < 1 is relevant, so that some fine tun- 
ing of the bare couplings fift is necessary to force jlf to 
flow into a fixed point of the RG. Because a finite limit 
/Et^ = lirnt_, 00 jlf means that we have found the true 
Fermi surface of the interacting system |J, we conclude 
that the detailed shape of the Fermi surface is sensitive 
to the numerical values of the bare couplings of the the- 
ory. Polchinski jll| pointed out that such a fine tuning of 
the bare couplings is unnatural, because physical effects 
which depend on the precise shape of the Fermi surface 
are not protected against small perturbations. 

To further elucidate the relation between the relevant 
couplings flf and the shape of the Fermi surface, let us 
explicitly derive from Eq. (|ll]) a self-consistency con- 
dition for the Fermi surface. Therefore it is useful to 
transform Eq. ( p"l| ) into an integral equation, 



jlf = e *-/o dT < 



dt'e-t'+ti' dr ^f[ 2) (n) 



(13) 



Suppose now that we have adjusted the bare couplings 
such that for t — > 00 the flowing couplings flf indeed 
approach finite fixed point values. Assuming that the 
associated anomalous dimensions rf^ are smaller than 
unity) 19 1 we conclude from Eq. (13) that the limit jj^ = 
lim^oo flf can only be finite if the initial values /i" are 
chosen such that 



Mo 



dte~ t+ $o dT ^tf\n) 



dte -t+tid T r,» I G t (Q')f( 4 )(Q ,Q';Q',Q ) 



(14) 

This is an implicit equation for /i", relating it to the 
values of the two-point vertex and the four-point vertex 
on the entire RG trajectory. Keeping in mind that the 



right-hand side of Eq. ( 
that according to Eq. ( | 



14) implicitly depends on jlf and 



E(k F ,iO) - E Co (k F ,zO) = &K/z$ 



(15) 



it is obvious that Eq. (|T^) can be regarded as an integral 
equation for the counterterm E(k F ,iO), the solution of 
which yields the true shape of the Fermi surface. At 
this point it is instructive to transform Eq. ( |l4| ) back 
to unrescaled variables, choosing for simplicity the initial 
conditions E£ o (k F ,?0) = and Z" = 1. Using the above 
definitions we find that Eq. (14) is equivalent with 



dk! du' 



) D 2tt iu3> - £ k ' +11- E 4k , {K') 



xTf (k F , iO, cr, k', iu', a'; k', iw', a', k F , iO, a) , (16) 

where £k' = |e k ' — £k' F |- Note that the right-hand side of 
Eq. ([h]) involves the flowing self-energy and four-point 
vertex at the scales £ = £k' which depend on the distance 
from the true Fermi surface. We emphasize that the exact 
integral equation (|l^) and the equivalent rescaled equa- 
tion ( |l4| ) fix the counterterm £(k F ,iO) from the require- 
ment that for t — > 00 all couplings approach finite fixed 
point values. 

If the system is a Luttinger liquid, then the rescaled 
version (14) of the self-consistency equation is more con- 
venient than the unrescaled version (|l(]), because for a 
Luttinger liquid the marginal part of the four-point ver- 
tex vertex without wavefunction renormalization fac- 
tors diverges for £ — > 0, while the rescaled four-point ver- 
tex defined in Eq.(|^) approaches for t — > 00 a finite 
limit. In this case the divergence of the unrescaled r^ 4 ' is 
canceled by the vanishing wavefunction renormalization 
factors at the Luttinger liquid fixed point p0| |. 

Given a solution E(k F ,iO) of Eq. (|l6|), we may calcu- 
late the compressibility of the system by substituting the 
result for E(k F ,iO) into Eq. (|]J) and solving for fc F (n), 
which implicitly depends on the chemical potential \i. 
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According to the Luttinger theorem Q the density n(/x) 
of the system is determined by the volume enclosed by 
the Fermi surface, 



n(M)=x; / ■ 



(Ik 



■0(Mn)-|k|) 



(2ir) D 

so that we obtain for the compressibility Xn 



n 2 Xr, 



dn x -> 



rf»5n D _x, M 9fcF(n) 
rip 



(2tt) d F 



9/i 



(17) 



(18) 



Note that the compressibility is a functional of the shape 
of the true Fermi surface of the many-body system, i.e. 
the compressibility characterizes the RG fixed point. Be- 
cause the rescaling of momenta and frequencies is crucial 
to obtain fixed points of the RG |L7| , the compressibility 
and other uniform susceptibilities arc not so easy to cal- 
culate within alternative versions of the RG which omit 
the rescaling 

From our point of view, the Fermi surface is a fixed 
'point manifold of the RG, so that it is meaningless to 
talk about the RG flow of the Fermi surface. However, 
it is possible to define a "flowing Fermi surface" k Ft via 



£k F , t - A* + £f oe -t(k.F )f ,iO) = 



(19) 



which by construction approaches the true Fermi surface 
for t — » oo. If we choose the initial conditions at t = 
(corresponding to £ = £o) such that S^ (k^o, iO) = 0, 
then kp\o is the Fermi surface of the non-interacting sys- 
tem at the same chemical potential as the interacting 
system. This corresponds to a non-interacting system at 
the density n (/ii) = I (27p e ( fc F,o(n) - |k|), which 
is in general different from the density n(fi) of the in- 
teracting system given in Eq. (|l7|). Note that in Fermi 
liquid theory one usually works at constant density. How- 
ever, as discussed by Nozieres ^l], for a self-consistent 
calculation of the Fermi surface it is more convenient to 
determine the counterterm at constant chemical poten- 
tial [A, and calculate the corresponding density afterwards 
p2[ . The practical advantages of such a procedure have 
been recognized previously in Ref. p3fl . Moreover, in 
the field-theoretical approach advanced by Ferraz the 
chemical potential is also a RG invariant. 

In the simplest approximation Eqs. (|l4|) and ( |l6| ) re- 
duce to the Hartree-Fock self-consistency equation for the 
counterterm E(ki?,iO). To see this, we approximate the 
flowing four-point vertex in Eq. ( |l6| ) as follows 

ri 4 ) (k F , iO, a, k', iv', a'; k', iw', a', k F ,iO, a) 

« r^ (k F , iO, cr, k' F , iO, a'; k' F , iO, a', kp, iO, a) 

= r< 4) (k f ,a;k F ,a') , (20) 

i.e. we project all momenta onto the Fermi surface, ignore 
the frequency dependence, and replace the flowing vertex 



by its fixed point value. Moreover, at this level of approx- 
imation we may also replace E^ k , (K') — > S(k' F , iO) on the 
right-hand side of Eq. (16). For £o - * °o we then obtain 
the Hartree-Fock self-consistency equation 

S(k Fj tO) = J2j 7^D r o 4) (kF,a;k' F ,a') 

x Gin-eu -£(k' F ,zO)) . (21) 

For a given interaction vertex T^\]s.p,a]k.' F ,a'), Eq. 
( pl| ) can be used to study possible Fermi surface insta- 
bilities such as the Pomeranchuk instability (24j at the 
Hartree-Fock level, which should always be the first step 
before more elaborate methods are used. A trivial gen- 
eralization of Eq. ( ^l| ) with a spin-dependent countert- 
erm £(kp, iO, a) leads to the self-consistent Hartree-Fock 
equation for spontaneous ferromagnetism, implying the 
usual Stoner instability for strong enough interactions in 
the spin-channel. More generally, if we allow for other 
types of symmetry breaking, within the same approxi- 
mations as above the exact RG fixed point equation for 
the two-point vertex can be reduced to the Hartree-Fock 
self-consistency equation for the corresponding order pa- 
rameter. Note that no approximation has been made 
in deriving Eqs. ( fl4] ) and Jl6|), so that these exact fixed 
point equations (or generalizations thereof for other types 
of symmetry breaking) can serve as a starting point for a 
systematic calculation of corrections to the Hartree-Fock 
approximation. 

In summary, in this work we have shown how the true 
Fermi surface of an interacting Fermi system can be de- 
fined self-consistently as a fixed point property of the RG. 
Our main result are the two equivalent integral equa- 
tions ( |l4| ) and (|l^), which determine the counterterm 
S(k F ,iO) necessary to calculate the Fermi surface from 
Eq. (0). In the simplest approximation, Eq. ( |l6| ) re- 
duces to the Hartree-Fock self-consistency condition for 
the counterterm. However, systematic improvements are 
possible. 

Very recently Dusuel and Dougot |2j| presented a de- 
tailed analysis of the Fermi surface deformations in quasi 
one-dimensional electronic systems, using perturbation 
theory and the RG method. They realized that some 
"slight modification" of the Wilson-Polchinski RG ap- 
proach is necessary in order to use this approach for a 
self-consistent calculation of the Fermi surface, but ad- 
mitted that a practical implementation of such a modi- 
fication remains to be attempted. We have shown here 
that the necessary modification of the functional RG used 
in Refs. [|l3|, [l4, 15, ll| is simply the usual rescaling step 
|1 , which is an essential part of the "orthodox" Wilsonian 
RG |l7j]. We conclude that the functional RG approach 
for interacting Fermi systems in the form advanced in 
Ref. H provides an elegant solution to the difficult 0] 
problem of self-consistently constructing the true Fermi 
surface. 
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